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1 Introduction 



Polynomials play key roles in many fields of the system theory [1], fundamental problems in auto- 
matic control, filter theory and network realization need to check some properties of polynomials, 
i and positivity of polynomials is an important one of such properties [2] . When checking positivity 

Q\ ' of polynomials using traditional methods for proving inequalities, complexities of algorithms are 

increasing rapidly as variable number increases [3]. Nowadays, Lu Yang [4] proposed a concise 
method to prove positivity of homogeneous polynomials (i.e., forms), that is difference sub- 
stitution (DS), or its varied form weighted difference substitution (WDS). This method 
demonstrates great efficiency. [5] further showed that if a form is indeed positive definite (PD) 
or not positive semidefinite (PSD), then these properties can be checked by finite steps of WDS. 
For integral forms, we estimate in this paper upper bounds of step numbers required in checking 
these properties, they only depend on the variable numbers, the degrees and the upper bounds of 
absolute values of coefficients of the forms. Therefore, we can also prove whether an integral form 
is PD through finite steps of WDS. 



2 Main Result 

We first introduce following definitions and notations according to [5]. 
Considering T„ <E R™ xn , where 



/ 1 - - \ 
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Denote by 8„ the set of all n! permutations of {1, 2, . . . , n}, for (k\k2 ■ ■ ■ k n ) G 0„, let -P^ /£,...&„) = 
(o-ij)nxn be the permutation corresponding (fcife ■ ■ ■ k n ), that is 

1, j = h 

o, i ± h ■ 

Let 

^(fcifca...fc n ) = P(k!k 2 ...k n )T n , 

and call it the WDS matrix determined by permutation (k\k,2---k„), denote by T n the set of all n\ 
such matrices. The variable substitution x = ^4(fe 1 fe 2 ...fe„)y corresponding {k\k% ■ ■ ■ k n ) is called a 
WDS, where x = (x\, x 2 , ■ ■ ■ , x n ) T , y = (j/i, J/2) ■ • • j 2/n) T , the following set of substitutions 

{x = A ± A 2 ■ ■ ■ A m y : A, eT n } 

is called the m-th WDS set. 

Suppose A = (dij) e R" x ™, we call it a normal matrix, and the corresponding substitution a 

n 

normal substitution if ^ a,y = 1, j = 1,2, . . . n. Thus WDS matrices are normal matrices and 
WDS substitutions are normal substitutions. 

Lemma 2.1. Let A = (aij) nX n = B1B2 ■ ■ ■ Bk, where Bi(i = 1, 2, . . . , k) are all normal matrices. 

n n 

Then A is a normal matrix, and for the substitution x = Ay, we have ^ Xi = j/,. 

i=l i=l 

Proof. Suppose £?i = (buj),B2 — {b 2 ij) are normal matrices, let C = B1B2, and denote by 
C = (cjj), then 

n nn n / n \ n 

i=i 1=1 fc=i fc=i \i=i / k=i 

Thus C is normal, and further we can prove A is normal by introduction. Moreover, we have 

n n n n / n \ n 

i=l i=l j=l j = l \i=l / j=l 



□ 



Let /(x) G R[xi, X2, ■ ■ ■ , x n ] be a form, we call 



WDS(/) = |J {f(A e x)} (2) 

the WDS set of /, 

WDSW(/)= |J ••• |J {/0V...A ei x)} (3) 

B m ee„ eiee„ 

the m-th WDS set of / for positive integer m, and set WDS (0) (/) = {/}. 

Denote by N the set of nonnegative integers, let a = (a±, a.2, ■ ■ . , a n ) £ N n , |a| = ai + «2 + 
• • • + a n . For a form of degree d 

f(x±,X2,...,X n )= S ' 1 x 2 2 • • ■ x n " , 

|a|=d 

if all coefficients c a are nonzero, we say / has complete monomials. 
Let 

R+ = {{xi,x 2 , ■ ■ ■ ,s„) T :xi>0,i = l,...,n}, 
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denote the (n — l)-dimcnsional simplex in R™ by 

A„ = |(xi,x 2 , . . . ,x n ) T ■.y^ j x i = 1, (x 1 ,x 2 , ■ ■ ■ ,x n ) T e R™ j . 

and let 

A„ = j(xi,x 2 , . . . ,x n ) T :^Txi<l, (xi,x 2 , ■ ■ ■ ,x n ) T e R™ j , 

Definition 2.1. Suppose D C R™, /(x) £ R[xi, x 2 , . . . , x n }, /(x) is PD in D if /(x) > /or any 
xefl\ {0}, and i* is PSD in D if /(x) > for anyxeD\ {0}. 

Obviously, we have, 

Lemma 2.2. .A form f S R[xi, x 2 , . . . , x n ] /ias i/*e sa?ne positivity in R™ , A„ and A„. 

Denote by Z the set of integers. We deduce the following result for integral forms. 

Theorem 2.1. Suppose f <E Z[xi, x 2 , . . . , x n ] is a form of degree d, and the absolute values of its 
coefficients do not exceed M , then we have 

1. f is PD in A„ , if and only if there exists m < C P (M, n, d), such that each form in WDS^(/) 
has complete monomials, and its coefficients are all positive; 

2. f is not PSD in A„ (i.e., the minimum of f in A„ is negative), if and only if there exists m < 
C nps {M,n,d), such that a form in WDS^"^(/) has complete monomials, and its coefficients 
are all negative. 

where 

"in ^2 d "M d " +1 n d " +1+d d(™ +1 ' d+,ld "(d + l)( ,l ~ 1 )("+ 2 )^ 



C p (M,n,d) = 



C nps (M,n,d) = 



Inn — ln(n — 1) 

In ^2 dn+1 M d " +1 n d " +1+d d ( - n+1 ' ,d+nd " (d + l)( ,l ^ 1 )( ?l + 2 )^ 
Inn — ln(n — 1) 



+ 2 (4) 



(5) 



Thus, we can completely determine positivity of / through checking positivity of coefficients 
of forms in WDS {C ^° {M < n ' d)) (/): 

1. If each form in WDS^ c '' lps ^ M ' n ' d '"(/) has complete monomials, and its coefficients are all 
positive, then / is PD in A„; 

2. If each form in WDS (c,lps(M '"' d)) (/) has a nonnegative coefficient, then / is PSD in A n ; 

3. If there exists a form in WDS'- c ''* I " , ^ Af '™' d ^(/) has complete monomials, and its coefficients 
are all negative, then / is not PSD in A„. 

3 Estimate for lower bounds of positive definite integral 
forms in the simplex 

7] gives estimate for lower bounds of positive definite integral polynomials in simplex, [s] improves 
the estimate. 

Lemma 3.1 ([8]). Suppose f £ Z[xi, x 2 , . . . , x n ] is positive definite in A n . If the degree of f is 
d, and absolute values of its coefficients do not exceed M , then 

mmf > (2M)- <i " +1 d-(" +1 ) d " +1 . (6) 

A„ 
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Indeed, the deduction in [8] has proved the following more general result. 

Lemma 3.2 ([8]). Suppose the minimum of f € Z[xi,X2, ■ ■ ■ ,x n ] in A„ is not zero. If the degree 
of f is d, and absolute values of its coefficients do not exceed M , then 



min/ 



> (2A/)- d " +1 d- [n+1)dn " 



(7) 



We have the following result for integral forms in A„. 

Lemma 3.3. Suppose the minimum of f G Z[xi,X2, ■ ■ ■ , x„] in A„ is not zero. If the degree of f 
is d, and absolute values of its coefficients do not exceed M , then 



min f 

A„ 



> Ci(M,n, d). 



(8) 



where Ci(M,n,d) = (2M)- d n 



—d n ^—d n —dj—nd™ 



i 



Proof. Let [x\ja, ■ ■ ■ , x n> o) be a minimal point of / in A„, then x^o > — ,1 < j < n, we can 

suppose x n o > — without loss of generality. Thus 

n 



min f 

A„ 



/(xi >0 , . . . ,X n>0 )| 

xi,a 



> 



(nx nfi ) 

Xl,0 



f 



X„-1,0 1 



nx n .o ' ' nx n .o ' n 



f 



Xn-1,0 1 



nx n fi nx n ,o n 



Lct 



g(xy, . . . ,x„_i) =n d f(x ll . . . ,x„_i, -) 



=f(nxi, . . . ,n.x„_i, 1), 

then its minimum is not zero in A„_i. Degree of p G Z[xi, . . . , x n _i] is d, and absolute values of 
its coefficients do not exceed n d M, so from Lemma 3.2, we have 



mm g 



> (2A/)- d "n~ d " +1 d-" d ". 



Since 



we have 



xi,o x„_i,q 



G A n _i, 



Therefore 



/ 



xi.o 



Zn-1,0 1 



nx n fi nx n fi n 



>n- d (2M)- dn n- dn d~ 
= {2M)- dn n- d " +1 ~ d d- ndn 

min/ > (2M)- d "n- dn+1 - d d- nd " 



□ 
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4 WDS and barycentric subdivision 

In the A n simplex coordinate system, considering a WDS 

x = T n y, (9) 

we can see that ai = (1, 0, 0) T is transformed to (1,0,..., 0) T , a.2 — (i, h, . ■ ■ , 0) T is transformed 
to (0, 1, . . . ,0) T , . . ., and a„ = (i, i, . . . , i) T is transformed to (0,0, ... , 1) T . Moreover, &k(k = 
1, 2, ■ • • , n) is the barycenter of the (k — l)-dimcnsional proper face containing aia2 • • • a& in A„. 
Since (9) is a normal substitution, from Lemma 2.1 we know, after transform (9), the corresponding 

n 

point for any (xi,x 2 , • • ■ ,x n ) T € A„ satisfies = 1, that is, coordinates after transforms are 

i=l 

also normal. So, aia2-a„ is a subsimplcx of A„ after the first barycentric subdivision, it 
corresponds a WDS matrix T n = (ai a2 • ■ ■ a n ). 

Similarly, other n! — 1 WDS matrices respectively correspond other n! — 1 subsimplcxes of 
A„ after the first barycentric subdivision. Thus, from geometrical views, a WDS corresponds a 
barycentric subdivision of A„. 

From Lemma 2.1 and the definition of WDS, we know that sequential WDS correspond se- 
quential barycentric subdivisions of A„. 

Denote by diam a the dimension of simplex a, i.e., maximal distance between vertexes of a. 
Comparing with the dimension of original simplex, dimensions of subsimplexes in barycentric 
subdivision decrease. That is 

Lemma 4.1 ([(>]). Let a be an n- dimensional simplex, a' is a subsimplex in the barycentric 
subdivision of a, then 

TL 

diam a' < diam a. (10) 

n + 1 

5 Proof of the main result 

Proof of the Theorem 2.1. We will prove two propositions respectively. 
(I) Sufficiency is obvious. Now we suppose / is positive definite in A„. 
Let 

f(x\i ■ ■ ■ j X n ) S Ci 1 ...i n Xi ■ ■ ■ X ™ , 

ii + . . 

choose an arbitrary m-th WDS for / 



xi = axyi + (ax + /3 12 )y 2 + ■ ■ . + (ax + /3x n )y n 
X2 = a 2 yi + (a 2 + ^22)2/2 + ■•• + (a 2 + /3 2 n)y n 

x n = a n yi + (a n + P n 2)y2 + • • - + (a n + /3 n 
where on, en + fyj > 0, i = 1, . . . , n, j = 2, . . . , n. Since a WDS is a normal substitution, 



(11) 



i=l 



From Lemma 4.1, we have 



Y>i<(^Y", i = 2 „, 



Further more. 



1^1 ^f 1 ^) , ' 1 »../ -2 (13) 
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Let u = yi + . . . + y n , then (11) can be written as 

xi = aiu + (3 12 y 2 + ■■■ + /3i n y v 
x 2 = a 2 u + (3 22 y 2 + ... + f3 2 „y n 



x n = a n u + (3 n2 y 2 + . . . + (3 nn y n 



So 



f {x\ , . . . , Xn) 

= f(aiu + (3\ 2 y 2 + . . . + f3 ln y n , . . . , a n u + (i n2 y 2 + . . . + p nn y n ) 

= E c ll ... in {aiu + p 12 y 2 + ... + f3 ln y n ) 11 ■■■{a n u + (i n2 y 2 + ... + (3 nn y n ) 1 

ii — d 

= E c n--i n a l ■ ■ ■ a n ud + E <fe-in( a l! • • • ) PrmWl ' ' ' Vn 
ii-\-...-\-i n — d ii-\-...-\-i n —d 

= /(«!,... ,a„)u d + E ^...^(ai, . . . (Aro)^ 1 • • 

ii + ...+i„=d 

= X! ( ■ I d ' ■ | /( ai ' <fc-i„(ai, . . .,0 ^I 1 • ■ 'Vn 



il+...+i„ 



E ^ii-inVT ■■■Vn, 



i\ + . . .+i n — d 

where 



and 



Y2 <t>ii-i n {ui, ■ ■ ■ , PnnWl ■■■Vl 
ii-\-...-\-i n —d 



E I E — j— -^l^VI ' ■ ■ Pln uni V2 2 ' ■ ■ Vn' 

l + -+in=- VilI + ...+i 1 »=i I Jll """ Jln " 



j 11 + ...+j rtl ^d 



E ■ r ^" 1 PI? ■ ■ ■ Pin V r ■ ■ ■ Vt 

, , ■ Jnl • Jnn • 

i-rii . . n n n 

e ^ e Tp^rn^nn^r- 

n+...+« n =d iu + .-- + ii„=n 1 1 P^9= 1 fc=l p=l g=2 

„jll + ---+3nl, Jl2+---+Jn2 ,,jln+- • •+3ni> 

T-rri . . n n n 

e ^ e rf^n^nn^ 

il + -- +i,i=d in + . . . = »! 11 P,« = 1 -'P'?' p=lg=2 

Jll+-'+Jjil^ d 

3nl + ■ ■ ■ + jnn = in 

E C?11 + ■ ■ ■ + JnlV- Sl j 12 + ...+jn2+S 2 . . „il„+...+i„„+S„ 
si! •••«„! 2/1 2/2 y " 

si+...+s„=jn+...+j„i 



P9 



(14) 



— —J(ai, . . . ,a n ) + ^...f (ai, . . . ,/?„„). (15) 

111 •••l n ' 
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from (12) and (13), we have 



t— rn . i n n n 
JP9- k=1 „_i o 



Cm + -..+Jnl)! 



<A/(d!) n+1 
=M(rf!) n+1 
<M(rf!) n+1 
<Md ( ™ +1)d 



p=l g =2 



n 



n 

n- 1 

n 
n — 1 

n 



There are 



d + n-1 
n — 1 



< (d+iy 



nonncgative integer tuples (ii , . . . , i n ) satisfying i\ + . . . + 
up by terms whose number does not exceed 

(d+l)(n-D(n+2) ) 

and the absolute value of each term does not exceed 

'n - 1 



d, so 4>i 1 ---i n (cui, . . . , /3 nn ) is summed 



Therefore 



|^ il ... <B (a 1> . .. ,0 nn ))\ <Md<" +1 >< , (d+l)<»- 1 ><"+ a > f^)' 



From Lemma 3.3, we have 



dl 



ii! ■■■i n \ 



/(ai,...,a„) > Ci(M,n,d). 



(16) 



(17) 



From (15), (16), (17), we know that in order that djk > 0, it suffices 



d(M,n,d) > Md {n+ ^ d {d+l) {n - 1){n+2) 



n — 1 



That is 



TO > 



In ^2 d "A/ d " +1 rt <i " +1+d d(" +1 ) (i+ncr (d + 1)(" _1 H™+ 2 )^ 



(18) 

inn — ln(n — 1) 

(II) Sufficiency is also obvious. Now we suppose the minimum of / in A„ is negative, and a 
minimal point is (ai, . . . , a n ). From Lemma 3.3, we have 

/(ai, . . . , a n )\ > Ci(M, n, d). 

Suppose (y n , . . . ,y ln ) T , (y 2 i, ■ ■ ■ ,2/2n) T G A„, (in, . . . ,x ln ) T ', (x 2 i, ■ ■ ■ 1 x 2n ) T are coordinates 
satisfying (11), from Lemma 2.1, we have (in, . . . , xi n ) T , (X21, ■ ■ ■ , X2n) T G A n . From the corre- 
spondence of WDS and barycentric subdivisions, we have 



/n-lV 

y/ (Xu - X 2 l) 2 + . . . + (x ln - X2„) 2 < ( — - — J 
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Let Sj = xij — X2j,j = 1, ■ ■ ■ ,n, then there exists 7 £ (0, 1), such that 



. . . ,Xl n ) - f(x 2 l, ■■■,X 2n 

+ ■ 

OXi 



+ ^"g^T ) /(^l + 7^1, ... , Z 2 n + 7<U 



< 



< 



< 



3=1 

n - 1 
n 
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O^-f {X21 +-ySi,... ,X 2n + jSn) 



J 

m n 



2J 2J \ c h-i n ( x 21 + "fSlY 1 ■ ■ ■ (X2j 
j=l ii+...+i„=d 
, n 

E E ^ 

3=1 ii+...-Mn=d 



(a?2n + 7<*t0* 



d _ 1 f d + n — 1\ /n — 1 



n 

n — 1 

71 



<nM2 d - 1 (d+l)™- 1 
Thus if m is sufficiently large, such that 

nM2 d - 1 (d+l) n - 1 

i.e., 



n — 1 



< -Ci(M,n,d), 



m > 



In (2 d " +d M d " +1 n d + rf + 1 d n < r (tf + l)"" 1 



m-th WDS (11) satisfies 



Inn — ln(n — 1) 



f(ai, . . .,a n ) > -Ci(M : n,d). 



From the deduction of (I), we can see that if 

In ^2 dn+1 ]\l dn+1 n dn+1+d d ( - n+1 ^ d+ndn (d + l)( n_1 )( n + 2 ) 



m > 



Inn — ln(n — 1) 



(19) 



(20) 
(21) 

(22) 



and m satisfies (20), there exists a form in WDS' ra '(/), it has complete monomials, and all 
coefficients are negative. Comparing the right hand sides of (20) and (22), the latter is larger, so 
it suffices (22). 

□ 
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